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Abstract 



We calculate one-loop renormalization factors of bilinear operators made of 
physical quark fields for domain-wall QCD. We find that finite parts of such 
renormalization factors have reasonable values at 1-loop except an overlap 
factor between the physical quark field and the zero mode in the theory. 
We point out that the 1-loop estimate of overall renormalization factors be- 
comes unreliable at the coupling where numerical simulations are currently 
performed, due to the presence of this overlap factor. We show that this prob- 
lem disappears if the mean-field improved perturbation theory is employed for 
renormalization factors. 
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I. INTRODUCTION 



The lack of chirally invariant fermion formulations is one of the most uncomfortable 
points theoretically and practically in lattice QCD. For example, in the Wilson fermion 
formulation, which is popularly used in numerical simulations, the chiral limit can be realized 
only by the fine tuning of bare mass parameter, which compensates the additive quantum 
correction to the quark mass. 

Recently the domain- wall fermion formulation which was originally proposed for 

lattice chiral gauge theories ||, has been employed in lattice QCD simulations Q and has 
shown its superiority over other formulations: there seems no need of the fine tuning to realize 
the chiral limit while there is no restriction to the number of flavors. In particular massless 
mode which presents at the tree level seems stable against the quantum correction. Since 
this property is not a trivial one, two of us have performed the perturbative calculation for 
the domain-wall quark propagator and have shown that the massless mode is indeed stable 
against the 1-loop correction ||. The wave function renormalization factor for the massless 
mode has also been evaluated there. 

In this paper we extend our previous perturbative calculation to an evaluation of renor- 
malization factors for quark bilinear operators at 1-loop in domain-wall QCD. These renor- 
malization factors are needed to convert matrix elements such as meson decay constants 
measured in lattice simulations to those in MS scheme, so that the values predicted on the 
lattice can be compared with the corresponding experimental values. In particular, in order 
to discuss the scaling behavior of such matrix elements, the correction by the renormaliza- 
tion factors becomes important, since 1-loop corrections are in general non- negligible at the 
gauge coupling used in simulations. Moreover since 1-loop corrections decrease only as the 
inverse of the logarithm of the lattice spacing, not a power of it, they can not be removed 
through the usual continuum extrapolation like linear or quadratic in the lattice spacing. 
Therefore the corrections by the renormalization factors to such matrix elements must be 
included for precise measurements in numerical simulations with the domain-wall QCD. 

We start from the Green's function with quark bilinear operators and physical quark 
fields. We show that the divergence in the Green's function can be renormalized into the 
quark propagator and bilinear operators. There appears no mixing between these operators 
and operators made of heavy unphysical fermions. The renormalization factor consists of 
two parts, one is the usual renormalization factor which has similar value as that of Wilson 
fermion, another is the factor of overlap between field on each two boundary wall and the 
zero mode. We notice that this overlap factor is renormalized additively and tends to be 
large. We present renormalization factors with mean-field improvement to suppress this 
large correction. 

This paper is organized as follows. In Sec. II we introduce the action of domain-wall 
QCD together with some notations used in this paper and explicitly give the Feynman rules 
relevant to our calculations. In Sec. Ill we evaluate the self-energy of the physical fermion, 
from which we extract the wave function renormalization factor of the physical fermion and 
the multiplicative mass renormalization factor. Although results for both wave function || 
and mass || renormalization factors have already been reported, we give them here in order 
to demonstrate our method of calculating renormalization factors in domain-wall QCD. In 
particular we pay attention to the extra factor in the wave function renormalization, which 
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arises as the overlap between the physical fermion field and the zero mode. We then calculate 
the renormalization factors for quark bilinear operators in Sec. IV. Numerical results of one- 
loop coefficients for various quantities are given in Sec. V as a function of M in both cases 
with and without the mean-field improvement. Our conclusion is given in Sec. VI. 

In this paper we set the lattice spacing a = 1 and take SU(N C ) gauge group with the 

- 1 

gauge coupling g and the second Casimir Cp 



2N, 



c 



II. ACTION AND FEYNMAN RULES 

The domain- wall fermion is a 4+1 dimensional Wilson fermion with a "mass term" which 
depends on the coordinate in the extra dimension. For the explicit form of the "mass term" 
we adopt the Shamir's one |]J in this paper and the quark action with the current quark 
mass m becomes 



JV 

Sew = 

n s=l 



\ Y {^{n) s {-r + 7^)t/ /t (n)^(n + /i) s + i[)(n) s (-r - 7m)^ (« ~ ~ I 1 )*) 

A ft 

+ - (^(n) s (l + 75 )V>(n) s+ i + ^(n).(l - 7 5 )^(n) a _i) + (M - 1 + 4r)^(n) s ^(n), 
+ mJ2 (?(n)ivi , +V'(n)i + ?(n)iP_V'(n)iyr) , (H-l) 

n 

where n is a 4 dimensional space-time coordinate and s is an extra fifth dimensional index, 
the Dirac "mass" M is a parameter of the theory which we set < M < 2 to realize the 
massless fermion at tree level, m is a physical quark mass, and the Wilson parameter is set 
to r = — 1. Our 7 matrix convention is as follows: 

7*=(£ 74= (J I), 75 = 71727374 = (J \ ) , (IL2) 

P± = ^4^, (H.3) 

^ = - [7m> >] • ( IL4 ) 

The gauge part of the domain- wall fermion is a standard 4 dimensional one plaquette action 
and we have no gauge interaction along the fifth dimension. We can interpret this fifth 
dimensional degrees of freedom as a flavor 0. It is important to notice that we have 
boundaries for the flavor space; 1 < s < N. 

The remarkable property of the domain-wall fermion is that there exists a massless 
fermion mode in the N — > oo limit at small momentum at m — 0. This massless fermion 
stays near the boundaries of the flavor space with the left and the right mode on the opposite 
side. At tree level the massless mode Xo is given explicitly in zero momentum limit as 



Xo = V 1 " (P+w'o-^s + P-<~ s ^s) , (IL5) 

where wo = 1 — M. Although the zero mode is stable against quantum correction f5j, the 
damping factor wq is modified because the Dirac mass M is renormalized additively. In 
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numerical simulations it is more convenient to use the interpolating "physical" quark field 
defined by the boundary fermions 



q(n) = P + ^( n )i + P-^{n) N , 
q(n) = ip(n) N P + + ip(n)iP— 



(II.6) 



We use the QCD operators constructed from this quark fields, since this field has been 
actually used in the previous simulations. Moreover, as seen later, we find that the renor- 
malization becomes simpler for the operators in terms of the physical quark field. 
Weak-coupling perturbation theory is developed by writing 

1 



U n ^ = exp(igaA^(n + -ft)). 



(11.7) 



Since our gauge part is same as that of the usual Wilson plaquette action, the gluon propa- 
gator can be written as 



G±(P) 



1 



4 sin p/2 



V - (1 



O: 



4sinp M /2 sinp^/2 
4 sin 2 p/2 



6, 



ab ■ 



(11.8) 



where sin 2 p/2 = X) A1 sin 2 p A1 /2 and we set a = 1 in our calculation. Quark-gluon vertices 
are also identical to those in the N flavor Wilson fermion. We employ the following two of 
those relevant for the one loop calculation: 



v i»(k,p)st = -2#T a {7 M cos^(-/^ +p^ - ir sin + Pn)}$st, 

V% v (k,p)st = ^ 2 i{T a ,T fc }{z 7At sin^(-A; M - rcos |(-^ +p v )}8 lu ,6 8t . 



(II.9) 

(n.io) 



Our momentum assignments for the vertices are given in Fig. [l|. 

The fermion propagator is given by inverting the domain-wall Dirac operator, whose 
formula in momentum space is 



f(p)st = J2^ sin P^st+(W + (p)+mM + ) st P + + {W-(p)+mM -) P-, (11.11) 



where the mass matrix is 



W + (p) 



s,t 



W-(p) s , t 

( 



f-W(p) 



\ 

f-W(p) 

1 



1 

-W(p) 



(11.12) 



-W(p) J 



■W(p) 



-W(p)J 



M~ 



W(p) = l-M-rJ^i 1 -cosp^). 



(11.13) 

(11.14) 
(11.15) 
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In numerical simulations one should take N large enough in order to have a massless quark in 
the spectrum at m = 0. This means that we can effectively neglect e~ aN terms with positive 
a. Therefore in our 1-loop calculation we will take N — > oo limit to avoid complications 
arising from the finite N. Although the N — > oo limit should be taken after the momentum 
integral in principle, it is easy to see that the limit can be taken before the integral. We 
may take the form of free fermion propagator in the N oo limit from the beginning. 
(See || for the derivation of the finite N fermion propagator.) Inverting ( [11.1 1| ) the fermion 
propagator for N —>■ oo becomes 



s F (p) st = (pipy 1 ) 



St 



= (-i^aupp + W- +mM + ) su G R (u,t)P + 
+ f-z7„ sinp, + W + + mM~) G L (u, t)P_, 

V / su 

where sum over the same index is taken implicitly. Gr/l is given by 
.4 



G R (s,t) 



[1 - m 2 ) (l - We' a ) e a{ ~ 2N+s+t) - (1 - m 2 ) (1 - We a ) e 



-2W sinh(a)m ( e *(- N + s -t) + e <-N-s+t) 



G L (s,t) = - 



F 



-a 



m 



(1 _ We ^ e a(-2N+s+t-2) _ (1 



+ Ae -a\s-t\^ 

1 - We~ a 



m 



cosh(a) 
A 
F 



-2W sinh(a)m ( e «(- N + s ~t) + e <*{-N-s+t) 

l + iy 2 + E M sin 2 p M 
2W 



+ Ae -a\s-t\ 



2M^sinh(a)' 
1 - e a W -m 2 [l-We 



(11.16) 

a{s+t) 

(11.17) 

e a(-s-t+2) 

(11.18) 

(11.19) 
(11.20) 
(11.21) 



Note that the argument p of factors a and W is suppressed in the above formula. 

Besides the fermion propagator given above, the propagator for the physical quark field, 
defined in ( [II .6|) , is also used in this paper. The explicit form is rather simple and is given 
by 



S q ip) = (q(-p)q(p)) 



-r^sinp^ + (1 - We' 



m 



(11.22) 



III. RENORMALIZATION FOR QUARK MASS AND WAVE FUNCTION 

In this section we calculate the self energy for the physical quark field, from which we 
derive renormalization factors for mass and wave function. The physical quark propagator 
( |11.22| ) at tree level becomes 

s g ( P ) = w \ (nil) 

if + (1 — Wo)m 
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in the continuum limit. In the following we will see that the full quark propagator with 
1-loop corrections takes the following form: 



full 



;i - wi)z w z 2 



+ wl)Z w Z- l m'' 



(III.2) 



where Z 2 and Z m are the quark wave function and the mass renormalization factors respec- 
tively. In Ref. || Z 2 is calculated, while Z m is given in Ref. ||. Z W) which is dropped in 
Ref. 0, is a renormalization factor for the overall factor (1 — Wq). 

One loop corrections to the quark two point function (q{—p)q{p)) are given by two 
diagrams in Fig. |2|. In the diagrams interaction vertices contain fermion field of a general 
flavor index, while physical quark field appears on the external lines. Tree-level propagators 
for the external lines are given by 



q(-p)ip(p,s 
1 



i (ij^smp^ -m (l - We~ a )) (e" a(7V - s) P + + e - a( ^ 1} P_ 



+ - 



m 



(ijn sinp^ — m (l — We 



& ~a / e - a ( S -l)p + + e -a(N-s)p_ 



m-p,s)q(p)) = i ( e -^)p_ + e~ a ^P + ) (i 1(i si 



smp^ — m 



(l - We' 



+- (e- Q(s - 1} P_ + e- a{N - s) P + ) e- a [m (ij^ sinp M - m (l - We 



-F 



(III.3) 



(III.4) 



In the one loop calculation we take external momenta and quark masses much smaller 
than the lattice cut-off, so that we can expand the propagator corresponding to external 
lines in terms of external momenta and masses. In diagrams in Fig. ^ it is enough for the 
renormalization to expand the external propagator to the next to leading order, 



l-wl 



+ 



w 



q(-p)ip(p,s i 

' \( (N-s) D . (s-1) D 

- ( W P+ + W Q P- 

■ (JUq J ill V 

(ip(-p,s)q(p)) 



Wo r+ 



Wo 



wl 



W P- 



1 



Wo)m 



(III.5) 



(III.6) 



Note that both terms at the leading and the next to leading orders in the external propagator 
give the leading order contribution in these diagrams since interaction vertices have the form 
a + bp^ due to the Wilson term. 

The external propagator has the characteristic form in the continuum limit that it con- 
sists of the quark propagator times some dumping factor. One loop correction to the quark 
propagator becomes 



(q{-p)q{p)) 1 



Wn 



Wo)m 



E q {p,m) 



WP, 



i4 + (l 



Wo)m 



(III.7) 



which is easily renormalized into the quark propagator. Here m) is given by multiplying 
the fermion self energy E(p, m) st by the factors in the external propagator, 
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E,(p,m) 



v ' v y 1 — Wn 



S(p,m) s t 

(III.8) 



Calculating the quark self-energy £(p, m) a t to one-loop order and making an expansion of 
the form 



£(/>.//»,, =-(()), + ^^^ + ^^m + 0(p 2 ,m 2 ,pm), 



9m 



(III.9) 



we find 



E(p, m) st = -g 2 C F [ij (/+P+ + I~P_) + W A 1 + P+ + Wf/L + m (m 1 +P+ + Mf P_ 



St ' 

111.10) 



where I^ 1 and are already calculated in our previous paper || for the massless fermion, 



r+A 



+ 



T 

2 st+ 16^' 
<i 4 Z 1 



(C+/-).,t ln(7r 2 )-l-lnA' 



(2vr) 4 4sin 2 1/2 
sin 2 L 



^4(4 sin 2 J/2) 2 



~ ^ (cos/ M (^-G fi + W + G L )(s,t) + sin 2 I^Gl + G R )(s,t)) 
+ W + G L )(s, t) + 2(2 cos 2 1 - 2 cos 2 t) 



-2T5 



st 



dH 



(27r) 4 4sin 2 //2 



- (C+/-).,* / 

E 



d 4 / 1 



(2vr) 4 (P) s 



7T 2 - I 2 ) 

cos 2 ±(W + /-G L/R )(s,t) 



fin.ir 



- sin 2 |(W-/+G fi/L )( S ,t) + ~sin 2 l^G L + G R ){s,t) 

Here T is the tadpole loop integral 

» 7r r/ 4 / 1 

0.15493. 



(111.12) 



tt (2vr) 4 4sin 



2 2 



(111.13) 



Note that the notation used in this paper is slightly different from that in Ref. ||. In 
particular the infrared divergence in the loop integral is regularized by a gluon mass A 
introduced in the gluon propagator in this paper, while the external fermion momentum 
is kept non-zero in Ref. || to avoid such an infrared divergence. The calculation becomes 
simpler in the former. 

It is also straightforward to calculate one loop correction to mass term Mf, 



(Mr 



+/- 



167T 5 



;P + /_) Slt (ln(7r 2 )-lnA< 



7 



+ 



dH 



(27r) 4 4sin 2 1/2 



E 



cos 



2 ^ (w +/ - oa 



|f + M + /- Gl/s )(M) 



— sm 



k ( W-'^ + M-'+G^ M) + i sin* J?^ + °£ 



dm 
dH 1 



dm dm 



where C+/_ and B+/_ are denned by 

{C + ) st = (1 - «; >r 2 , (C_) s , = (1 - Wo 2 H 2JV — ' 



(B + ) st = (1 - 



2\2„„7V-s+i-l 



(S_) s , = (1 - wlfw* 



2\2„„s-l+N-t 



(s,t) 
(111.14) 



(111.15) 
(111.16) 



Substituting these into ( |III.8|) and neglecting higher order terms of 0(p 2 ,m 2 ,pm) we 
obtain 



1 -u9 2 Cf 



— m- 



1 — Wq 167T 2 

g2CF (-A log A 



(-logAV - Si 



2u>n 



16tt 2 



2^2 v 



where 



,T 



S 1 = -log7r 2 + l + (|--/ d )-167r 2 , 

E 2 = -41og7r 2 -a m -167r 2 , 
S 3 = (2T - D) • 16vr 2 , 



(111.17) 



(111.18) 

(111.19) 
(111.20) 



+ E 



/ (S { 32^ S H + ^ + 2 C ° S ^ - ^* 

(w - + (£ cos 2 Z„/2 - 2 cos 2 IJ2)G L + $>in 2 Z„/2)G R 

V V 

(111.21) 



sin Z M 
2(4 sin 2 Z/2) 2 



(Z 2 ) 2 
ci 4 / 



UJ 



vr 2 - / 2 ) 
1 



(27r) 4 4sin 2 //2 



E 



cos 2 |( Wo - W)G R - sin 2 |( Wo - W)G R 



+ -sin 2 / M (G L + G R ) 
rf 4 / 1 



(111.22) 



o. 



(27r) 4 4sin 2 //2(l- U ; e- a ) 



E 



cos 



2.2'^ —2a , • 2 ; 

+ r sm — e + r sm /, 



2" ' ■ "~ M l-e°iy_ 
The propagators and Gr are defined by 



, 2 Z„ 1 - e~ a V^ 
2 1 - e a W 

d 4 Z 1 



(2tt) 4 (^ 



7T 2 -/ 2 ). 



(111.23) 
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Gi 



G 



R 



2Wsmha 
2Wsmha 



sinh ao — sinh a 



e a -W 



1 



2wq sinh «o(cosh «o — cosh a) e~ a — W(e a — wq) 2 

sinh ao — sinh a 1 
2w sinh a (cosh a — cosh a) (e a — w ) 2 



(111.24) 
(111.25) 



with e a ° = Wq. It should be noted that there is no additive mass correction in Eq. ([111.171) . 
Although the fermion self energy E(p, m) s t has an additive mass correction, it vanishes in 



S 9 because of the relation, 



which has been shown in the previous paper [H [] . 

The full quark propagator at 1-loop level now becomes 



Sq(p) 



full 



i$ + (1 — Wo)m 



1 + E q (p,m) 



l-wl 



+ (1 - u>o) m - (1 - w o)^q(p, m) 



+ (1 — Wg)m 
+ 0(g A ), 



and finally we obtain 



full 



;i-™ 2 )z w z 2 



zjf + (1 - wl)Z w Zj-m 



with 



z 2 


= 1 + 


Z w 


= 1 - 




= 1 + 



9 2 C F 
16vr 2 
2w g 2 C F 



(logAV + Ej 
E 3 = l 



<? 2 C F 



1 — 167T 2 " 167T 

16tt 2 



2 



-31ogAV -£ 2 + Ei) ■ 



From this relation we can see that E 3 is an additive renormalization to Wq, 

9 2 C F ^ \ 2 



1 - w ) Z w = 1 - \Wq + 



16tt 2 



If we employ the mean-field (tadpole) improvement ||, we can write 

(i - (^nzfii-^ 



# + (1 - W F ) 2 )ZMF( M ZMF)-l m ' 



(111.26) 



(111.27) 



(111.28) 

(111.29) 
(111.30) 
(111.31) 



(111.32) 



(111.33) 



1 At finite N we need a few extra terms for the renormalization such as an additive mass counter 
term. However these terms always contain the factor e~ aN and can be suppressed in the large TV" 
limit before taking the continuum limit. 
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where u = 1 - g 2 C F T/2, w^ F = w - 4{u - 1), 

logA 2 a 2 + £ a - 16tt 2 ^ 



^•MF 


= 1 + 


9 2 C F / 

167T 2 V 


log A 2 


7 MF 


= 1 - 




>C F 




16tt 2 




= 1 + 


9 2 C F / 

167T 2 I 


-31o 



16vr 2 • 2T = 1 + 



-31ogA 2 a 2 -S 2 + S 1 -167r- 



9 2 Cf mf 

167T 2 10 ' 



(111.34) 
(111.35) 
(111.36) 



Of course equalities that uZ 2 = Z 2 MF , (1 - w%)Z w = (1 - (w^ F ) 2 )Z^ F and Z m = uZ™ F hold 
up to 0(g 4 ) terms. 

In the continuum we employ the MS scheme with naive dimensional regularization. The 
one-loop self-energy in the continuum has the same form as ( |IH.17|) with, however, the 
replacements, 



log(Aa) 2 - MA///) 2 , 

£^£^ = 1/2, 

V _^ vMS _ o 
Zj2 — ► ^2 — — Z ' 

w -> 0. 



(111.37) 
(111.38) 
(111.39) 
(111.40) 



Let us define the quark wave function renormalization factor needed for converting the 
lattice field to the continuum field in the MS scheme by 



g MS = (i-»or i/2 ^ 1/2 ^H? ,at - 

To one-loop order we finally find that 



Z 2 {\id) = 1 + -zr^C F - log(yua) 2 + z 2 



where 



and Z w is given in (|IH.30| ). The quark mass renormalization factor defined by 

m MS (/i) = (l - wt) Z w Z m (jj,a)m 



is given by 



Z m (iia) = 1 + -pr^C F -31og(/za) 2 + z n 



;m.4i) 



(111.42) 



(111.43) 



(111.44) 



(111.45) 



with 



MS - s 2 ) 



(Ef s - Ex] 



(111.46) 



10 



IV. RENORMALIZATION FOR QUARK BILINEAR OPERATORS 



We consider quark bilinear operators in the following form 

°t(x) = q(x)Tq(x) } T = 1, 75, 7^, 7^75, 0>. 



(IV.l) 



We calculate the one loop correction to the Green's function (Or(x)q(y)q(z)) for massless 
quark with external momenta p = p' = 0, whose diagram is shown in Fig. || As in the 
previous section we notice that the external line is essentially written in terms of the quark 
propagator times dumping factors, ( 111.5 ) and ( |111.6| ). Making use of this fact the one loop 
full Green's function for small external momentum becomes 



{(qTq) ■ qq) m 
where Tr is given by 

T r = 1 + 



1 - tug) Z W Z 2 _ T r (1 - wl) Z W Z 2 



(IV.2) 



9 2 C, 
16tt 2 



log X 2 a 2 + V T 



(IV.3) 



with a T dependent constant h 2 (T) = 4(A), 4(V), 16(P), 16(S), 0(T). The finite renormal- 
ization factor Vp is given by 



It 



h 2 (T) 



log 7T 2 + I r ■ 16vr 2 



(IV.4) 



dH 
h 2 (T) 



2 o I srn2 1 
r sin — h r 

2 e~ a -W 



+ 



.4, 



Wf\ (e Q - w ) Z 4sin 2 | 



dH 



(2vr) 4 (/ 



1X1 



(IV.5) 



where 



9 / • 2 1 

COS 2 sm ' 



r = s, p 

cos 2 I - 2 cos 2 l -f) (sin 2 I - 2 sin 2 

cos 2 I - 2 cos 2 l f - 2 cos 2 |) (sin 2 / - 2 sin 2 l a - 2 sin 2 ^) T = T a 



r — v^,, v4 Q 



(IV.6) 



sin 2 / = X^si 1 * 2 ^ an d cos 2 1/2 = cos 2 / /t /2. The renormalization factor of quark bilinear 
operator is written as 



1 - w 2 )Z w Z 2 T r = (1 - w )^r 



lat 



16tt 2 



log A 2 a 2 + V T + £1 



(IV.7) 
(IV.8) 



It is seen from ( [IV. 6| ) that Za = Zy and Zp = Z$, which is expected to hold if the chiral 
symmetry exists and is indeed so in the continuum QCD. 

In the continuum, the on-shell vertex function to one-loop order is given in the MS 
scheme by 
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t ms 1 , 9 2 Cf fh(T) 



1 + 



16tt 2 



MS 



(IV.9) 



with Vp MS = -1/2(A), -1/2(1/), 2(P), 2(5), 0(T) for the ant i- commuting definition of 75. 

Combining the above results and including self-energy corrections, the relation between 
the continuum operator in the MS scheme and the lattice operator is given by 



O r MS M = (1 



w, 



where 



Zr(na) = 1 + 



9 2 Ct 
16vr 2 



0) X Z 

II! 

Z r (fia) (qTq) , 
^-l) log(H 2 + zr 



and 



with 



.MS 



lat 



V MS 



^ r MS + 1/2 = 0(A), 0(F), 5/2(P), 5/2(5), 1/2(T), 



(IV.10) 



(IV.11) 



(IV.12) 



(IV.13) 
(IV.14) 



Here it should be mentioned that z m = —z s and consequently Z m (/ia) = Z s (^a) 1 , since 
a m — Is ( S 2 = — Vs )• To show this we use the relation 



[e- a - W) {e a -W) 1 + W 2 -2W cosh a 



sin / 



{e- a - W) 2 



[e~ a - wy 



-wy 



(IV.15) 



Note that the chiral Ward-Takahashi identities, if exact, implies the three equalities z$ 
zp = —z m and zy = z^, which is satisfied without fine tuning in the domain- wall QCD. 



V. NUMERICAL RESULTS 

To calculate momentum integrals appeared in the previous two sections numerically, two 
independent methods are employed, in order to check formulae and programs. In the first 
method the momentum integration is performed by a mode sum for a periodic box of a size 
L 4 after transforming the momentum variable through — g M — sin q^. We employ the size 
L = 64 for integrals. In the second method the momentum integration is carried out by the 
Monte Carlo routine VEGAS, using 20 samples of 1000000 points each. We found that the 
both result agrees very well. 

We first present basic quantities S 1; S 3 , Vs(— Vp), Vy(= Va), and Vr in table | for 
several values of M. Here S 2 is omitted since S 2 = —Vs holds exactly. The finite part of 
quark renormalization factors z w and Z2, are given in table 0. Since z m = —zs we omit 
it here again. We also give mean-field improved values where M should be replaced by 
M = M + A{u — 1). In table [ITJ we give our main result of this paper, zs,p, Zy t A and zt, for 
several values of M. In these tables errors estimated from difference between two methods 
of the numerical integration are in the (right) next to the last digit of each value. 
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VI. CONCLUSION 



In this paper we have calculated renormalization factors for quark bilinear operators 
in domain-wall QCD at 1-loop of the perturbation theory. We find that z$ = zp = —z m 
and zy = za, which suggest that the chiral Ward-Takahashi identity holds exactly without 
fine tuning at this order of the perturbation theory. Our numerical values for zx with 
X = 2, m, S, P, V, A, T are compared with those for the Wilson fermion action with and 
without mean-field improvement and clover fermion action with c$w — 1> given in table [IV]. 
No data for clover fermion action with mean-field improvement is given here since c$w itself 
is modified through the mean- field improvement. Finite parts zx are rather large without 
mean-field improvement in all three fermion formulations, but they becomes smaller after 
the mean-field improvement is performed for the domain-wall fermion and Wilson fermion. 

Peculiar feature of the domain-wall QCD in its renormalization is an appearance of the 
overlap factor (1 — Wq)Z w for the physical quark field. The problem here is that the 1-loop 
correction z w becomes very large unless |1 — M\ is small. For example, at M — 1.7, which 
is the value used in the quenched simulation at /3 = 6.0 ||, from table pi] z w = 137.03 
without mean-field improvement. This 1-loop correction is huge, so we can not trust the 
perturbation theory. If mean-field improvement is made with u = P 1 / 4 where P is the 
expectation value of the plaquette and is 0.59374 at (3 — 6.0, z w = 0.1893 at M ~ 1.2. This 
value is reasonably small. This example suggests that one have to employ the mean-field 
improvement for perturbative renormalization factors in domain-wall QCD to convert the 
lattice result to the one in the MS scheme. 
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TABLES 



TABLE I. Value of Si, S 3 and V T . 



M 


Si 


s 3 


V s ,p 


Vv,A 


V T 


0.050 


13.25 


51.22 


3.30 


4.833 


5.344 


0.100 


13.16 


51.05 


3.82 


4.835 


5.173 


0.150 


13.08 


50.89 


4.24 


4.836 


5.035 


0.200 


13.01 


50.74 


4.60 


4.837 


4.915 


0.250 


12.94 


50.61 


4.93 


4.838 


4.807 


0.300 


12.88 


50.49 


5.24 


4.840 


4.707 


0.350 


12.83 


50.37 


5.53 


4.84 


4.613 


0.400 


12.77 


50.27 


5.80 


4.84 


4.524 


0.450 


12.72 


50.17 


6.06 


4.85 


4.439 


0.500 


12.68 


50.07 


6.31 


4.85 


4.357 


0.550 


12.63 


49.98 


6.56 


4.85 


4.277 


0.600 


12.59 


49.90 


6.80 


4.85 


4.199 


0.650 


12.56 


49.83 


7.04 


4.85 


4.122 


0.700 


12.52 


49.76 


7.28 


4.85 


4.045 


0.750 


12.49 


49.69 


7.52 


4.86 


3.970 


0.800 


12.46 


49.63 


7.75 


4.86 


3.894 


0.850 


12.44 


49.58 


7.99 


4.86 


3.819 


0.900 


12.41 


49.53 


8.23 


4.86 


3.743 


0.950 


12.39 


49.49 


8.47 


4.87 


3.666 


1.000 


12.38 


49.46 


8.71 


4.87 


3.588 


1.050 


12.36 


49.43 


8.96 


4.87 


3.509 


1.100 


12.35 


49.41 


9.22 


4.88 


3.429 


1.150 


12.34 


49.39 


9.48 


4.88 


3.346 


1.200 


12.34 


49.39 


9.75 


4.88 


3.261 


1.250 


12.34 


49.39 


10.02 


4.89 


3.173 


1.300 


12.35 


49.40 


10.31 


4.89 


3.082 


1.350 


12.36 


49.42 


10.62 


4.89 


2.987 


1.400 


12.37 


49.45 


10.93 


4.90 


2.888 


1.450 


12.39 


49.49 


11.27 


4.90 


2.783 


1.500 


12.42 


49.54 


11.62 


4.91 


2.673 


1.550 


12.45 


49.61 


12.00 


4.92 


2.555 


1.600 


12.49 


49.69 


12.40 


4.92 


2.428 


1.650 


12.54 


49.80 


12.84 


4.93 


2.292 


1.700 


12.60 


49.92 


13.31 


4.94 


2.142 


1.750 


12.68 


50.07 


13.84 


4.94 


1.978 


1.800 


12.77 


50.25 


14.42 


4.95 


1.793 


1.850 


12.87 


50.46 


15.09 


4.96 


1.582 


1.900 


13.00 


50.72 


15.88 


4.97 


1.334 


1.950 


13.15 


51.03 


16.84 


4.98 


1.027 
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TABLE II. Value of zi and z w . Mean- field improved values are also listed. 



M 


Z2 


Zyj 


Z2 


MF 


Zw 


0.050 


-12.75 


-998.23 


-0.52 




-44.70 


0.100 


-12.66 


-483.60 


-0.43 




-20.05 


0.150 


-12.58 


-311.75 


-0.35 




-11.99 


0.200 


-12.51 


-225.53 


-0.28 




-8.057 


0.250 


-12.44 


-173.52 


-0.21 




-5.759 


0.300 


-12.38 


-138.59 


-0.15 




-4.272 


0.350 


-12.33 


-113.39 


-0.09 




-3.244 


0.400 


-12.27 


-94.248 


-0.04 




-2.502 


0.450 


-12.22 


-79.114 


0.01 




-1.946 


0.500 


-12.18 


-66.762 


0.06 




-1.521 


0.550 


-12.13 


-56.408 


0.10 




-1.188 


0.600 


-12.09 


-47.526 


0.14 




-0.9253 


0.650 


-12.06 


-39.748 


0.18 




-0.7147 


0.700 


-12.02 


-32.807 


0.21 




-0.5446 


0.750 


-11.99 


-26.503 


0.24 




-0.4062 


0.800 


-11.96 


-20.681 


0.27 




-0.2929 


0.850 


-11.94 


-15.217 


0.30 




-0.1995 


0.900 


-11.91 


-10.007 


0.32 




-0.1218 


0.950 


-11.89 


-4.9617 


0.34 




-0.05631 


1.000 


-11.88 


0.0 


0.36 




0.0 


1.050 


-11.86 


4.9553 


0.37 




0.04993 


1.100 


-11.85 


9.9813 


0.38 




0.09620 


1.150 


-11.84 


15.159 


0.39 




0.1416 


1.200 


-11.84 


20.577 


0.39 




0.1893 


1.250 


-11.84 


26.339 


0.39 




0.2428 


1.300 


-11.85 


32.569 


0.39 




0.3066 


1.350 


-11.86 


39.420 


0.38 




0.3863 


1.400 


-11.87 


47.091 


0.36 




0.4894 


1.450 


-11.89 


55.846 


0.34 




0.6264 


1.500 


-11.92 


66.053 


0.32 




0.8119 


1.550 


-11.95 


78.235 


0.28 




1.068 


1.600 


-11.99 


93.173 


0.24 




1.427 


1 650 


-1 9 04 


1 1 9 09 


1Q 




1 Q45 


1.700 


-12.10 


137.03 


0.13 




2.711 


1.750 


-12.18 


171.66 


0.05 




3.895 


1.800 


-12.27 


223.31 


-0.03 




5.840 


1.850 


-12.37 


309.12 


-0.14 




9.357 


1.900 


-12.50 


480.47 


-0.26 




16.91 


1.950 


-12.65 


994.48 


-0.42 




40.95 
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TABLE III. Value of zr- Mean-field improved values are also listed. 



M 


ZS,P = -Z m 


ZV,A 


ZT 


ZS,P = -Z m 


MF 

ZV,A 


ZT 


0.050 


-14.05 


-18.083 


-18.09 


-1.82 


-5.8497 


-5.86 


0.100 


-14.48 


-17.996 


-17.83 


-2.25 


-5.7627 


-5.60 


0.150 


-14.82 


-17.918 


-17.62 


-2.59 


-5.6853 


-5.38 


0.200 


-15.11 


-17.848 


-17.43 


-2.88 


-5.6149 


-5.19 


0.250 


-15.38 


-17.783 


-17.25 


-3.15 


-5.5503 


-5.02 


0.300 


-15.62 


-17.723 


-17.09 


-3.39 


-5.4906 


-4.86 


0.350 


-15.85 


-17.668 


-16.94 


-3.62 


-5.4352 


-4.71 


0.400 


-16.07 


-17.617 


-16.80 


-3.84 


-5.3838 


-4.57 


0.450 


-16.28 


-17.569 


-16.66 


-4.05 


-5.3359 


-4.43 


0.500 


-16.49 


-17.524 


-16.53 


-4.26 


-5.2913 


-4.30 


0.550 


-16.70 


-17.483 


-16.41 


-4.46 


-5.2500 


-4.18 


0.600 


-16.90 


-17.444 


-16.29 


-4.67 


-5.2117 


-4.06 


0.650 


-17.10 


-17.409 


-16.18 


-4.87 


-5.1764 


-3.95 


0.700 


-17.30 


-17.377 


-16.07 


-5.07 


-5.1440 


-3.84 


0.750 


-17.51 


-17.347 


-15.96 


-5.27 


-5.1145 


-3.73 


0.800 


-17.71 


-17.321 


-15.86 


-5.48 


-5.0879 


-3.62 


0.850 


-17.92 


-17.297 


-15.75 


-5.69 


-5.0643 


-3.52 


0.900 


-18.14 


-17.277 


-15.66 


-5.91 


-5.0437 


-3.42 


0.950 


-18.36 


-17.259 


-15.56 


-6.13 


-5.0262 


-3.33 


1.000 


-18.59 


-17.245 


-15.46 


-6.35 


-5.0119 


-3.23 


1.050 


-18.82 


-17.234 


-15.37 


-6.59 


-5.0012 


-3.14 


1.100 


-19.07 


-17.226 


-15.28 


-6.84 


-4.9935 


-3.05 


1.150 


-19.32 


-17.223 


-15.19 


-7.09 


-4.9898 


-2.96 


1.200 


-19.59 


-17.223 


-15.10 


-7.35 


-4.9900 


-2.87 


1.250 


-19.87 


-17.227 


-15.01 


-7.63 


-4.9944 


-2.78 


1.300 


-20.16 


-17.236 


-14.93 


-7.93 


-5.0034 


-2.70 


1.350 


-20.47 


-17.250 


-14.84 


-8.24 


-5.0174 


-2.61 


1.400 


-20.80 


-17.270 


-14.76 


-8.57 


-5.0368 


-2.53 


1.450 


-21.16 


-17.295 


-14.67 


-8.92 


-5.0621 


-2.44 


1.500 


-21.54 


-17.327 


-14.59 


-9.30 


-5.0941 


-2.36 


1.550 


-21.95 


-17.366 


-14.51 


-9.71 


-5.1335 


-2.27 


1.600 


-22.39 


-17.414 


-14.42 


-10.16 


-5.181 


-2.19 


1 650 


-99 88 


-1 7 471 


-14 34 


-1 65 


-5 938 


-9 10 


1.700 


-23.42 


-17.539 


-14.25 


-11.18 


-5.307 


-2.01 


1.750 


-24.01 


-17.620 


-14.16 


-11.78 


-5.388 


-1.92 


1.800 


-24.69 


-17.716 


-14.06 


-12.46 


-5.484 


-1.83 


1.850 


-25.46 


-17.831 


-13.95 


-13.23 


-5.598 


-1.72 


1.900 


-26.37 


-17.968 


-13.83 


-14.14 


-5.735 


-1.60 


1.950 


-27.49 


-18.135 


-13.68 


-15.26 


-5.902 


-1.45 
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TABLE IV. Value of zx for Wilson and clover fermion actions. 



Wilson Clover(c 5M / = 1) Wilson(MF) 

~Z2 -12.852 -9.206 -0.619 

z m 12.953 19.311 0.719 

z s -12.953 12.661 -0.72 

z P -22.596 9.602 -10.363 

z v -20.618 16.657 -8.385 

z A -15.797 18.193 -3.564 

z T -17.018 20.824 -4.785 
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FIGURES 




FIG. 1. Momentum assignment of the fermion-gluon vertices. 




FIG. 2. One loop correction to the quark propagator. 



qTq 




FIG. 3. One loop correction to the quark bilinear operator. 
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